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Abstract. There is an interesting proposal that the long-range Coulomb interaction 
in suspended graphene can generate a dynamical gap, which leads to a semimetal- 
insulator phase transition. We revisit this problem by solving the self-consistent Dyson- 
Schwinger equations of wave function renormalization and fermion gap. In order to 
satisfy the Ward identity, a suitable vertex function is introduced. The impacts of 
singular velocity renormalization and dynamical screening on gap generation are both 
included in this formalism, and prove to be very important. We obtain a critical 
interaction strength, 3.2 < a c < 3.3, which is larger than the physical value a = 2.16 
for suspended graphene. It therefore turns out that suspended graphene is a semimetal, 
rather than insulator, at zero temperature. 
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1. Introduction 

After the successful fabrication of monolayer graphene in laboratory p], tremendous 
experimental and theoretical efforts have been devoted to exploring its novel and 
fascinating properties [21 |3j HI |5] . Without any interactions, the low energy quasiparticle 
excitations of graphene are massless Dirac fermions with linear dispersion [6]. The 
long-range Coulomb interaction between Dirac fermions is poorly screened due to the 
vanishing of density of states vanishes at Dirac points, and is therefore anticipated to 
be responsible for many unusual behaviors of graphene [21 HI [3]. Interestingly, it can 
lead to singular renormalization of fermion velocity [Tl El El HQ] , which is able to induce 
unconventional properties in several observable quantities, such as specific heat [Tlj [T2]. 
compressibility [121 [13], and electrical conductivity [HI [151 EH E]- R is remarkable 
that the predicted singular velocity renormalization has already been observed in recent 
experiments [18] . 

If the long-range Coulomb interaction is sufficiently strong, a finite fermion gap 
may be dynamically generated through the formation of excitonic particle-hole pairs 

[191 Eol EH [221 Eal [241 EHl [261 EZ1 EHl EHJ Eol EU [321 [331 E4l EHl EE1 ET1 EH! [391 Sol ST]. 

The dynamical gap fundamentally changes the ground state of graphene [12] ■ As a 
consequence, there will be a quantum phase transition from semimetal to excitonic 
insulator. This gap generating mechanism is a concrete realization of the non- 
perturbative phenomenon of dynamical chiral symmetry breaking that has been 
extensively investigated in particle physics for five decades (131 0U 05J EH H7]. From 
a technological point of view, a finite gap would make graphene more promising as a 
candidate material for manipulating novel electronic devices [121 EE]. For these reasons, 
the dynamical gap generation and the resultant semimetal-insulator transition have 
attracted intense theoretical interest in recent years. Many analytical and computational 
tools, including Dyson-Schwinger(DS) gap equation [191 EQl EU [221 E31 EH ESI ES ET1 
[23129], Bethe-Salpeter equation [301 [31], renormalization group (RG) [321 [331 [31] , an d 
large scale Monte Carlo simulation [35l [36l [371 E3 ESI SO], have been applied to study 
this problem. A critical interaction strength a c and a critical fermion flavor N c are 
found in almost all these investigations: a dynamical gap is opened only when a > a c 
and N < N c . If one fixes the physical flavor N = 2, the semimetal-insulator transition 
is turned solely by the parameter a, with a c defining the quantum critical point. We 
list some existing values of a c for N = 2 in Table I. With a few exceptions, a c obtained 
in most calculations is smaller than the physical value a = 2.16 of graphene suspended 
in vacuum. Therefore, it is suggested by many that the suspended graphene should be 
an insulator at zero temperature. 

However, there is little experimental evidence for the predicted insulating ground 
state. Actually, a recent experiment put an upper limit, as small as O.lmeV, on the 
possible gap in suspended graphene [18]. It is known that a finite gap is observed in 
graphene which is placed on specific substrate [49] or has finite-size configurations [50J. 
Nevertheless, this gap appears to be induced by the very particular environments, and 
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can not be regarded as a true dynamical gap generated by chiral vacuum condensation. 
Generally speaking, there may be several reasons for the lack of clear experimental 
evidence for dynamical gap. For instance, the dynamical gap can be strongly suppressed 
by thermal fluctuation pH EH ED E2], doping [22], and/or disorders [221 123, which 
makes it technically hard to measure the dynamical gap in realistic experiments. 
Another possibility is that the Coulomb interaction is simply not strong enough to 
open a dynamical gap even in the clean and zero-temperature limits. 

Table 1. Existing predictions for the critical interaction strength at N = 2. To 
make a comparison between the results obtained in different references, here we define 
a = e 2 Ivf£. 



Reference 




Reference 


Ctc 


m m 


2.33 


[29] 


1.02 


m 


1.1 


m 


1.62 


m 


1.2 


[32] 


0.833 


m 


1.13 


[33] 


2.5 


m 


0.92 


m 


0.99 


m\ 


8.1 


[35] 


1.11(6) 


1281 


1.79 


1391 


1.66 



Motivated by the recent theoretical and experimental progress, we revisit this 
problem in order to specify the genuine ground state of suspended graphene. Here we 
only consider clean suspended graphene at zero doping and zero temperature. The 
gap generation will be examined by means of DS equation, which is known to be 
a very powerful tool of analyzing dynamical gap generation in a number of strongly 
interacting models [131 0U 05] , such as Nambu- Jona-Lasinio model [131 HI] , quantum 
chromodynamics (QCD) [H1H5], three-dimensional quantum electrodynamics (QED 3 ) 
pf6"l Wl\ l5Tj . and graphene. 

In Ref. [19], Khveshchenko applied the DS equation to examine the possibility 
of dynamical gap generation in graphene. In order to simplify the very complicated 
nonlinear gap equation, a number of approximations are introduced. First of all, the 
energy dependence of the polarization function is neglected, which is usually called 
instantaneous approximation in literature [IHl [20l Ell [221 [231 Ell ESI [2TI EH] - In addition, 
the renormalizations of wave function and fermion velocity are both ignored. In the 
subsequent years, the existence of dynamical gap generation in suspended graphene was 
rechecked after improving some of these approximations [221 E31 EH EH] , and it became 
clear that these neglected effects can significantly increase or decrease ct c . Unfortunately, 
due to the formal complexity of DS equation(s), it often happens that one specific 
approximation is improved at the cost of introducing another. In Refs. [231 ES], the 
effect of fermion velocity renormalization on gap generation is investigated and showed to 
increase ct c . However, the dynamical part of polarization function is not well addressed 
in their calculations. In Ref. [22] , a full dynamical polarization function is utilized, but 
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the velocity renormalization is not included. Ref. [24J performs a detailed analysis 
of the influence of dynamical polarization function without including the velocity 
renormalization and the energy-dependence of dynamical gap. Another potentially 
important effect is the strong fermion damping caused by Coulomb interaction. It 
is predicted that the long-ranged Coulomb interaction, even though being too weak to 
open a gap, produces marginal Fermi liquid behavior [91 [191 E2l [53]. The strong fermion 
damping may affect dynamical gap generation, but is rarely considered in the existing 
literature. Since a precise value of a c is crucial to answer the question regarding the 
true ground state of graphene, it is necessary to calculate a c after going beyond all these 
approximations. 

We will show that all the aforementioned effects can be incorporated self- 
consistently by constructing a set of DS equations of wave function renormalization 
A),iG°o ? p) and dynamical fermion gap m(p Q ,p). To satisfy the Ward identity, we also 
include the vertex correction to the fermion self-energy. The dynamical fermion gap 
and fermion velocity renormalization can be obtained simultaneously after solving these 
equations. Our numerical calculations lead to a critical coupling 3.2 < a c < 3.3 for 
physical flavor N = 2. Apparently, a c is larger than a = 2.16, so the ground state 
of suspended graphene seems to be a semimetal, rather than an insulator. Our a c 
is quite different from those obtained in previous DS equation studies, which implies 
that an appropriate treatment of velocity renormalization and dynamical screening is 
very important. Moreover, from A 0t i(p ,p), we recover the singular renormalization of 
fermion velocity. 

The rest of the paper is organized as follows. In Sec. 2, we construct and numerically 
solve the self-consistent DS equations. A critical interaction strength a c is obtained from 
the solutions. In Sec. 3, we address the feedback effect of fermion velocity renormalization 
on the effective Coulomb interaction. We briefly summarize our results and discuss the 
implications in Sec. 4. 

2. Dyson-Schwinger equations for gap and wave functions 

The Hamiltonian for interacting Dirac fermions is 



where the density operator Pi(r) = ipi(r)'yoijji(r) and g = 2-7re 2 /vpe. The Dirac fermions 
have totally eight indices: two sublattices, two spins, two valleys. As usual, we adopt 
a four-component spinor field ip to describe the Dirac fermions [191 120] , and define its 
conjugate field as ip = V> 7o- Now the physical flavor of fermions is N = 2. The 4x4 
7-matrices, 70,1,2, satisfy the Clifford algebra JTSJ [2U]. It is easy to check that the total 
Hamiltonian possesses a continuous chiral symmetry, ip — > e 1015 ^, where the matrix 
75 anticommutes with 70,1,2- If the massless Dirac fermions acquire a finite mass gap 
due to excitonic condensation, (ipip) 7^ 0, then this continuous chiral symmetry will be 



- Jr 4Jz7T . , / r r / r I" 





(1) 
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dynamically broken. 

It is convenient to characterize the interaction strength by an effective fine structure 
constant, a = e 2 /vp£- We need to consider the strong-coupling regime since the 
dynamical gap can only be opened by strong interaction. When a is close to or larger 
than unity, the conventional perturbation expansion in terms of a breaks down and one 
has to use 1/N as the expansion parameter. 




Figure 1. (a) The DS equation for effective Coulomb interaction, (b) The DS equation 
for fermion propagator. 



The free propagator for massless Dirac fermion is 

G (po,k) = , . (2) 

«7oPo - vf! • P 

Due to the Coulomb interaction, this propagator will receive self-energy corrections and 
become 

G(po, p) = —— ; v x r, (3) 

iA (p , p)7oPo - v F Ai(p , p)7 • p - m(p , p) 

where Ao,i(po, p) are the wave function renormalization and m(po, p) is the fermion gap 
function. Here we would like to make a comparison with the gap generation problem 
of QED 3 . In QED 3 , the explicit Lorentz invariance ensures that A (p ,p) = Ai(p ,p) 
[4"6| Wf\ . On the contrary, the Lorentz invariance is explicitly broken in the present 
system, so A Q (p ,p) ^ Ai(p ,p) and the fermion velocity is renormalized. The 
renormalized velocity is given by Ai(p , p)/A (p , p). The function A (p ,p) itself is 
also important since it determines the fermion damping rate caused by the Coulomb 
interaction. The dynamical fermion gap is described by m(p Q , p), which acquires a finite 
value when excitonic particle-hole bound states are formed due to sufficiently strong 
interaction. In previous DS equation analysis, a common approximation is to simply 
assume that, A (p ,p) = Ai(p ,p) = 1, which drastically simplify the complicated 
DS equations. However, the effects of velocity renormalization and fermion damping 
can not be well addressed by doing so. In order to include these effects in a self- 
consistent manner, it is better to analyze the coupled equations of Aq{j>q, p), Ai(p , p), 
and m(p , p). 

According to the Feynman diagram shown in Fig. 1(b), the relationship between the 
free and renormalized fermion propagators is formally determined by the following DS 
equation, 

G^ipo, p) = Go 1 (po, p) + J ^ (Ir) 2 1 " ^ ' P ' k ° 7 k ' )G ^ ' k ^° V (P° ~ fc o, P - k), (4) 
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where r (po, P; k , k) is the full vertex function and V(p Q — fc , p — k) is the effective 
Coulomb interaction function. Regarding the Coulomb potential as the time component 
of a gauge potential, one can obtain a Ward identity |9J that connects the fermion 
propagator and the vertex function, 

r (po,p;po,p) = — ^ • (5) 

idpa 

Apparently, To = 70 only when Aq(po, p) = 1. As we are willing to examine the effects 
of wave function renormalization on dynamical gap generation, the vertex function Tq 
can not be simply replaced by the bare vertex 70. In principle, one could build an 
equation of T (po,p; k Q ,k) and couple it self-consistently to Eq.(4). However, this will 
make the problem intractable. A more practical way is to assume a proper Ansatz. 
Choosing a suitable vertex function is a highly nontrivial problem, and have been 
investigated extensively in the contexts of various quantum field theories [15]. The 
research experience that has been accumulated in QED 3 [51] is especially helpful since 
this model is very similar in structure to our present model. There are several frequently 
used Ansatze for the vertex function [51]. Among these Ansatze, for computational 
convenience we choose the following one, 

T (po, p; fco, k) = / (p , p; k , k) 70 = - [A (p , p) + A (k , k)] 7 . (6) 

Other possible Ansatze of r (po, P; ko, k) can be analyzed by the same procedure. 

Substituting the full fermion propagator G(po, p) into the above DS equation and 
then taking trace on both sides, we can obtain the equation for dynamical gap m(po, p). 
In order to derive the equations of A (p ,p) and Ai(p ,p), we multiply both sides of 
Eq.fll]) by 70P0 and 7 • p respectively, and then take trace on both sides. After these 
manipulations, we finally arrive in a set of self-consistently coupled integral equations, 

1 f +oo dk f d 2 k [A (p ,p) + A (k ,k)]Mko,k)k 
A o{Po, P) — J- + - — / — 



2p Q J 2tt J (2n) 2 A 2 (k ,k)k 2 + v 2 F A 2 (k ,k^\ 2 + rn 2 (k ,k) 

x [V(p +k ,p-k)- V(p -ko,p- k)] , (7) 

1 [ + °°dk f d 2 k [A (p ,p) + A (k ,k)]A l (k ,k)(p-k) 

M,P) + 2\p\ 2 J 2n J (2n) 2 A 2 (k ,k)k 2 + v 2 F A 2 (k M^\ 2 + rn 2 (k ,k) 

x [V(p + k ,p-k) + V(p -k ,p- k)] , (8) 

. 1 f +oo dk f d 2 k [A (po,p) + A (A; ,k)]m(A;o,k) 

m{p , P) - 



2j 2tt J (2n) 2 A 2 (k ,k)k 2 + v 2 F A 2 (k Mk\ 2 + m 2 (k ,k) 

x [V( Po + fc , p - k) + V(p - fco, P - k)] . (9) 

During the derivation of these equations, we have used the following symmetries 

A),iG°o, p) = Au(-po, p), (10) 

m(po,p) =m(-p ,p), (11) 
which originate from the particle-hole symmetry of undoped graphene. 
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We now consider the effective interaction function V(po — k ,p — k). The bare 
Coulomb interaction is 

T , . . 2nav p , s 

(q) = "TqT' 

2 

with effective fine structure constant, a = —. The value of dielectric constant e depends 
on the substrate on which the graphene is placed. For graphene suspended in vacuum, 
the corresponding value is roughly a = 2.16; for graphene placed on Si0 2 substrate, 
a = 0.79 [12] ■ Besides the static screening due to substrate, the Coulomb interaction is 
also dynamically screened by the collective particle-hole excitations. After including the 
dynamical screening, as displayed in Fig. ([Ik), the effective interaction function becomes 

V(qo, q) = —17 y, (13) 

K) (q) + n(g ,q) 

The dynamical polarization function 11(g) is defined as 

/d 3 k 
-^^Tr [70^0(^0, q)7oG (A;o + go, k + q)] (14) 

to the leading order, which amounts to random phase approximation (RPA). It is 
straightforward to obtain 



N |q 



2 



n(g ,q) = -^=^==. (15) 
8 V9o + v f\<1\ 

Apparently, although the bare Coulomb interaction Vo(q) is independent of energy, 
the effective interaction function V(qo, q) depend on both energy q and momentum q 
due to dynamical screening. This V(q ,q) should be substituted to the self-consistent 
equations (7-9). 

Before doing numerical calculations, it is interesting to compare the present problem 
with QED 3 . In QED 3 , the inverse of bare photon propagator is ex g 2 , with g being 
three-dimensional energy-momentum. The corresponding polarization is known to be 
oc |g|, which is larger than q 2 at low energies [46J. Therefore, the effective photon 
propagator is dominated by the polarization and the bare term is relatively unimportant. 
In the present problem, the inverse of bare Coulomb interaction is Vg -1 (q) oc |q|. The 
polarization II (go, q) is more complex than that in QED3, but is evidently linear in 
momentum |q|. Therefore, the bare term Vo _1 (q) and the polarization n(g , q) appearing 
in the denominator of effective interaction V(qo, q) are equally important, which is 
different from QED3. 

After time-consuming but straightforward numerical calculations, we obtain the 
relationship between m(0, 0) and a, shown in Fig.Q. It is difficult to get a precise a c 
because the necessary computational time becomes increasingly long as a is approaching 
its critical value. Therefore, here we can only give a narrow range of a c . For physical 
fermion flavor N = 2, the critical strength is 3.2 < a c < 3.3. This critical value is 
obviously larger than a = 2.16, so it turns out that the zero-temperature ground state 
of suspended graphene is a semimetal in the clean limit. This result is consistent with 
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Figure 2. Relationship between dynamical gap m(0, 0) and a. Apparently, m(0, 0) 
decreases as a decreases. m(0, 0) becomes zero once a decreases below certain critical 
value a c . It is technically hard to get a precise a c . Numerical calculations show that 
m(0, 0) remains finite at a = 3.3 but vanishes at a = 3.2, so the critical value must be 
3.2 < a c < 3.3. 




for the graphene placed on SiC>2 substrate with a = 0.79. 



the fact that so far no experimental evidence for insulating ground state of suspended 
graphene has been reported. 

Our critical value a c is very different from the values presented in earlier 
publications. Compared with the previous gap equation computations, the key 
improvements in our analysis are that we have included the wave function 
renormalization Aq^(pq,p) self-consistently and maintained the energy- dependence of 
all the functions of A 0) i(p , p), m(p , p), and n(g , q)- We also have adopted an Ansatz 
for the vertex function in order not to violate the Ward identity. Our results indicate 
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that fermion velocity renormalization, fermion damping, and dynamical screening of 
Coulomb interaction all play important roles in determining the critical interaction 
strength a c . 

It is now useful to make a more concrete comparison between our a c and those 
presented in some recent literature. For example, our a c is much larger than the a c 
obtained by Liu et al. [22] and Gamayun et al. [23], who considered a fully dynamical 
polarization but ignored wave functions and velocity renormalization. Sabio et. al. 
[25] studied dynamical gap generation by means of a variational method that naturally 
includes velocity renormalization, and found a critical value a c = 8.1, which is much 
greater than our a c . Such a large quantitative difference is presumably owing to the 
instantaneous approximation of Coulomb interaction adopted in their analysis. Once 
the instantaneous approximation is assumed, the effective Coulomb interaction no longer 
depends on energy, and V(po + ko, p — k) will be identically equal to V(po — k , p — k). 
As shown in Eq.(7), the time component of wave function becomes A (p , p) = 1. Now 
the fermion damping effect is automatically neglected and the velocity renormalization 
is solely determined by Ai(p ,p). Moreover, A (p ,p) = 1 implies that vertex function 
is simply Tq = 70. After these simplifications, our DS equations become similar to those 
presented in Ref. [25]. By solving these simplified DS equations, we find no evidence for 
dynamical gap even for infinite coupling a — » 00. This indicates that the instantaneous 
approximation misses very important fluctuation effects. However, despite the large 
difference in the magnitude of a c , our conclusion do agree with Sabio et al. that no 
dynamical gap is generated in suspended graphene. 




Figure 4. (a) Renormalized velocity Vp(po,p) in the po —> limit for different values 
of a; (b) Renormalized velocity Vp(po,p) in the po — > 00 limit for different values of 
o. 

Although the Coulomb interaction in suspended graphene is too weak to generate 
a dynamical gap, it still leads to unusual properties: strong fermion damping and 
singular velocity renormalization. Usually, the fermion damping rate and the fermion 
velocity renormalization are calculated separately. Our DS equation approach includes 
the mutual influence of these two quantities self-consistently, so the damping rate and 
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velocity renormalization can be simultaneously obtained from the solutions of A (p , p) 
and Ai(po,p). The energy and momentum dependence of A (p ,p) and Ai(p ,p) are 
shown in Fig. ((3]), with (a-b) for suspended graphene with a = 2.16 and (d-e) for 
graphene placed on Si0 2 substrate with a = 0.79. It is clear that both A (p ,p) 
and A\(po, p) increase as interaction strength a is increasing. A Q (p ,p) decreases with 
growing energy and momentum. Different from A (p ,p), Ai(p ,p) is a decreasing 
function of momentum |p| but an increasing function of energy p . When |p| decreases 
gradually Ai(p , p) keeps increasing monotonically; whenp - > or p — > oo, A±(po, |p|) 
saturates to a finite value Ai(0,p) or Ai(oo,p), respectively. 

The fermion damping is determined by the time component of wave function 
renormalization, A (p ,p). For any given energy p , A (p ,p) increases as momentum 
|p| decreases in the region |p| > po, but saturates to a finite value A (p ,0) as |p| 
decreases in the region |p| < p . For given momentum |p|, when p decreases gradually, 
A (p , p) increases with decreasing p in the region p > |p|, but it saturates to a finite 
value v4o(0,p) with decreasing po i n the region po < \p\. In the low energy regime, 
A (p ,p) can be approximated as A (p ,p) ~ 1 + F(A/ a/Po + P 2 )' where F(x) is an 
increasing function of x. In principle, the fermion damping rate can be obtained from 
A (p ,p) by making analytic continuation, ip — » u> + i5. Unfortunately, the function 
F(A/ a/pq + p 2 ) is formally quite complicated and can not be written in terms of a 
simple analytical formula. Our numerical results suggest that F(A/ \/pl + p 2 ) grows a 
little more slowly than log(A/ \Jp$ + P 2 ) as \lv\ + P 2 decreases continuously, hence the 
damping rate does not display an exact linear dependence on energy. By ignoring the 
singular velocity renormalization, some pervious perturbative [9j [52] and self-consistent 
[19l 153] calculations predict a marginal Fermi liquid behavior, namely T(u) ~ oj. Such 
linear-in-energy behavior disappears once the singular velocity renormalization is self- 
consistently considered. However, the renormalization factor Z extracted from our 
numerical results is found to vanish, Z = 0, so the suspended graphene does not have 
well-defined quasiparticles. In this sense, our results are qualitatively consistent with 
the previous conclusions [EJ [52j [53] . 

The fermion velocity renormalization is related to many properties of graphene 
[3 El [H HOl [HI H21 UHl [HI UHl [16l [17]. Such effect is usually studied by RG method 
[3 EJ El ITU] . It can be naturally obtained by solving the self-consistent DS equations. 
The ratio between renormalized and bare fermion velocities is defined as 

Vf(Po, P) = A(po,p) ^ 
v F A (p , p) ' 

We plot Vp(po, p) for a = 2.16 and a = 0.79 in (c) and (f) of Fig. ([3]), respectively. The 
renormalized fermion velocity Vp(po, p) increases singularly as |p| — > 0, which reproduces 
the previous results of RG analysis [TJ 13 [10] . We also obtain the energy dependence 
of renormalzied velocity Vp(po, p), which is not well captured in RG analysis. For given 
momentum, Vp(p , p) is an increasing function of energy. When p — )• or p — >• oo, 
v f{Pqi p) ceases to change, but is saturated to Vp(0, p) or Vp(oo, p), as can be seen from 
Fig. (3c) and Fig.(3f) as well as Fig. (4). 
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Figure 5. (a), (b), (c) and (d) are m(p , \p\),A (po, |p|), A 1 (p ,\p\) and v§(p ,\p\) 
for the graphene with a = 4.0. 

It is also interesting to consider the gapped phase with a > a c , though no known 
graphene material manifests such a large a. When a = 4.0, a dynamical gap is generated 
by the Coulomb interaction. The corresponding m(jo ,p), A 0t i(p ,p), and Vp(p ,q) are 
shown in Fig. ([5]). The gap m(p ,p) decreases with growing momentum but increases 
with growing energy. These behaviors are very different from those in QED3 where the 
dynamical gap is a decreasing function of both energy and momentum. This difference 
reflects the different energy- dependence of Coulomb interaction and gauge interaction. 
From Eq. (jl~3]) . we know that the effective Coulomb interaction decreases with growing 
momentum but increases with growing energy. However, the gauge interaction in QED3 
always decreases as energy or momentum grows, which is in accordance with the fact 
that QED 3 exhibits asymptotic freedom [4*6j W7\. 

In the low-energy regime, \/p\ + |p| 2 < m(0, 0), A(p 0} p) saturates to a finite value 
and can be approximate as A (p , p) ~ 1 + F(A/ a/W) IpI 2 m2 (0> 0)) with F(x) being 
certain increasing function of x. After analytic continuation, ipo — > 00 + i5, the fermion 
damping rate is found to vanish when \oj\ < a/|p| 2 + m 2 (0, 0). It implies that the 
fermions can not be excited below the scale of m(0,0). Unlike the case of semimetal 
phase, v§(po, p) does not keep increasing as momentum decreases. For any given energy, 
v f(Po^p) saturates to Vp(p ,0) below the scale |p| ~ m(0, 0), so the singular velocity 
renormalization is suppressed by the dynamical gap generation in the insulating phase, 
as shown in Fig. (4). These results are consistent with those obtained by perturbative 
calculations [54] and by functional RG calculations [55] . 
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3. Effects of velocity renormalization on polarization 

In the calculations presented above, we have used the RPA expression of the dynamical 
polarization function II (g , q), in which the fermion velocity is a constant. However, this 
function may receive sizeable feedback effects from the velocity renormalization. Since 
the effective Coulomb interaction plays a crucial role in determining a c , it deserves 
to examine these effects. In a fully self-consistent analysis, one should consider the 
following defination 

/d 3 k 
j^Tt [T G(k , k) lo G(k + q ,k + q)] , (17) 

where G(k ,k) is the full fermion propagator given by Eq.(3) and T is the vertex 
function. By doing so, the polarization n(g , q) couples self-consistently to the equations 
of Ao,i(po)P) an d m(po,p). Unfortunately, this would make the computational time 
unacceptably long since the integrations over energy and momenta must be performed 
separately in the present problem. To simplify the calculations, we assume that the 
RPA expression of n(g , q) provides a reliable description of the dynamical screening 
effect, which then allows us to replace the constant velocity Vp appearing in n(g , q) 
by the renormalized velocity Vp(p ,p). This strategy is analogous to that employed in 
Ref. 1 231. After this substitution, the effective interaction becomes 



V(q ,q) 
It is larger than 



1 



|q| , N |qr 

2naVF 8 v^o+KC^q)) 2 ^! 2 



|q| [_ N |q| 



! " (19) 



2 ™ VF 8 v^o+^lll 2 

since Vp(qo, q) > vp, so the velocity renormalization tends to enhance the effective 
interaction strength. In order to examine the extend to which this feedback effect 
changes a c , we re-solve the coupled gap equations after substituting Eq.(18) into Eqs.(7 
-9). We find the fermion gap is always zero when a < 2.9. It implies that such feedback 
effect does not change a c substantially. In particular, a c obtained after considering this 
effect is still larger than the physical value a = 2.16 in suspended graphene. 



4. Summary and discussions 

In conclusion, we have presented a detailed analysis of dynamic gap generation due 
to Coulomb interaction in graphene using the coupled DS equations of wave function 
renormalization and fermion gap. After including the effects of fermion velocity 
renormalization and dynamical screening of Coulomb interaction, we obtain a critical 
interaction strength, 3.2 < a c < 3.3. It is apparently greater than the physical strength 
a = 2.16 in suspended graphene. Therefore, the Coulomb interaction in suspended 
graphene is too weak to generate a dynamical fermion gap, and the semimetal ground 
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state of suspended graphene is robust, which is consistent with the fact that so far no 
insulating phase has been observed in experiments. 

Although our calculations suggest that the Coulomb interaction in suspended 
graphene is not strong enough to open a dynamical gap, the possibility of dynamical 
fermion gap generation can not be entirely precluded. Indeed, there do exist several 
possible mechanisms for the dynamical generation of fermion gap. For instance, 
the Dirac fermions can acquire a tiny bare gap for some reasons, such as Kekule 
distortion [56] and spin-orbit coupling [57]. When this happens, a large dynamical 
fermion gap can be generated by a relatively weak Coulomb interaction [26], which then 
produces properties analogous to the remarkable phenomena of QCD [26J. Moreover, 
the additional on-site repulsive interaction may help to generate a dynamical gap even 
if the Coulomb interaction itself is not sufficiently strong j22j [24j [27]. Finally, it is well 
known that an external magnetic field perpendicular to the graphene plane can lead to 
dynamical gap generation and insulating phase even at infinitesimal coupling [58| [20] . 
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